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A model 'remarkable' fin equation is singled out from a class of nonlinear (1 + l)-dimensional 
fin equations. For this equation a number of exact solutions are constructed by means of using 
both classical Lie algorithm and different modern techniques (functional separation of variables, 
' generalized conditional symmetries, hidden symmetries etc). 
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1 Introduction 



Mathematical models of conductivity and diffusion processes are traditional objects for investi- 
gations with symmetry methods and related to them approaches jlU| . Since these models are 
often formulated in terms of nonlinear differential equations which are, as rule, non-integrable 
^ ' and cannot be linearized, symmetry methods are important for construction of their exact so- 

lutions. It is indicative that modern development of group analysis of differential equations 
, was begun from the study of a class of (1 -|- l)-dimensional nonlinear diffusion equations jl7j . 

I Afterward a wide range of diffusion equations was investigated within the symmetry framework 

^ ■ (see e.g. [21 13 El El El US] ) • At the same time, even 'simplified' (1 -|- l)-dimensional non-linear 

, diffusion models are fraught with a great many of 'symmetry mysteries' which remain to be 

^ \ solved. 

C^' Recently the class of nonlinear fin equations of the general form 

■ ut = {D{u)u^)x + h{x)u, (1) 



was investigated with the symmetry point of view in a number of papers (31^11213 125 • Here u is 
treated as the dimensionless temperature, t and x the dimensionless time and space variables, D 
the thermal conductivity, h = —N'^f{x), N the fin parameter and / the heat transfer coefficient. 
[ (See e.g. for references on physical meaning and applications of equations (0).) 

Note that certain partial equations from class were studied formerly. For example, the 
condition = corresponds to the linear case of which was completely investigated with 
the Lie symmetry point of view long time ago |14l I18j . The problem of group classification 
for the class of nonlinear one-dimensional diffusion equations (the degenerate case h = 0) was 
first solved by Ovsiannikov jl7l I18j . The class of diffusion-reaction equations classified by 
Dorodnitsyn jH El contains the equations of form with h = const. Group classification of 
the subclass where the thermal conductivity is a power function of the temperature was carried 
out in 123] • Large sets of exact solutions were constructed for the above equations and collected 
e.g. in UniEIl. 

In contrast to |31 El 1201 ) the study in |211 was concentrated on rigorous and exhaustive group 
classification of the whole class and on construction of exact solutions for truly nonlinear 
and 'variable-coefficient' equations from this class. To find exact solutions, both classical Lie 
reduction and different modern approaches were applied. Although some interesting exact so- 
lutions were constructed, almost all of them are either stationary or scale- invariant. Therefore, 
the problem on more complicated exact solutions of equations with D^hx ^ remains still 
open even for partial values of the parameter-functions. 
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In this letter we single out a model 'remarkable' fin equation with the fixed values of the 
parameter-functions D = u^^/'^ and h = i.e. 

ut = {\u\~^/'^u-x)^ + x~'^u, (2) 

and investigate it in detail. A number of its exact solutions are constructed with a variety of 
symmetry techniques (Lie reduction, non-linear separation of variables, generalized conditional 
symmetries, hidden symmetries etc) in a closed form. 

2 Lie invariance 

Equation is ordinary, from the Lie symmetry point of view, in class (^. The maximal Lie 
invariance group of © is Ai = D = tdt + xdx — ^udu), i.e. its Lie symmetry group Gi 
consists of the transformations 

i=e^H + 6o, x = e^^x, u = e~^^^u, 

where 6q and 5i are arbitrary constants. At the same time, equation (jSJ has remarkable proper- 
ties connected with different kind of non-Lie symmetries that allows us to construct a number 
of its exact solutions. Sources of this singularity should be investigated additionally. 

Due to possibility of changing sign of u and due to physical sense of the equation, we can 
assume u positive and omit modular in the expression |n|~'^/^. 

Instead of equation 0, we can investigate the equivalent equation 

Vt = OTxx - - (3) 

6 2 X 

where v = u"'^!"^ and, therefore, v is positive. In the variables (t, x, v) the operator dt has the 
same form and D = tdt + '^dx + vd„. The inverse transformation from ((31) to © is u = v~'^l'^ . 

3 List of exact solutions 

For usability we collect all constructed solutions of Q together and then discuss methods of 
finding them, including both Lie and non-Lie techniques. The adduced solutions are inequivalent 
with respect to the group Gi and can be extended to parametric sets of solutions of equation Q 
with transformations from this group: 



/ 9 \-2/3 

\) u=\ -E^x^ Sex^ ~ 4^ ) ' ^ ^ ■L"'^' ^' ^J' '^1' 



9 x-2/3 



2) u=\ -X 

' \2 t ^ 

3) = {^^ - "ix^ tan 2t - 

4) w = ^-x^ 3x2 ^^^Yi 2t - 

5) u= ^-x^ + 3x2 pQ^j^ 2t - 

, 1 1 

6) VV'-V'2--arcsin(2^-l) = ±--FCo, := , < V < 1, 

2 XX 

7) VV^+V^-7tM2V' + 1 + VV^+V^) = ±- + C'o, V:= , ^/;<-lorV>0. 

2 XX 



2/3 
-2/3 
-2/3 



2 



Solutions l)-5) should be considered only for the values of {t, x) where the corresponding bases 
of power —2/3 is positive. 

Note that solutions of © can be transformed to solutions of equations which are pointwise 
equivalent to equation ((21). For example, the transformation i = t, x = , u = x'^u links (jSJ 
with the equation 

More generally, the above equation and (UJ belong to the class of variable-coefficient diffusion- 
reaction equations which were investigated in |23]. Therefore, equation ((2)) can be extended by 
transformations from the corresponding equivalence group to a subclass of variable-coefficient 
diffusion-reaction equations with simultaneous extension of the above solutions. 



4 Lie reductions 

Let us discuss ways of finding the above solutions starting with the classical Lie method. Ai is 
a non-Abelian two-dimensional Lie algebra. A complete list of inequivalent subalgebras of Ai is 
exhausted by the one-dimensional subalgebras {dt) and {D) and the algebra Ai itself. 

With Ai we construct the ansatz u = Lpx~'^^^ . In view of positivity of n, f should be 
also positive. Since Ai has a single functionally independent invariant, there are no invariant 
independent variables in the ansatz and (p \s a, 0-ary function i.e. a constant. Therefore, the 
ansatz reduces equation ^ to the algebraic equation Aip"'^/"^ + 9 sign x = with respect to Lp, 
which has a solution only on the negative semiaxis x < 0. As a result, we obtain solution 1) 
with e = 0. It is interesting in the sense that all other solutions l)-5) are modifications of it 
with additional terms. 

The solutions being invariant with respect to the subalgebra {dt) are nothing but stationary 
solutions. The corresponding anzatz u = p{uj), where to = x, gives the reduced ODE 

+ c^" V = (4) 

which is integrable in quadratures. Equation (^ is connected with equations 6.101 and 6.205 
of Integrability of (Q can be explained in the framework of symmetry approach. The Lie 
invariance algebra of (Q is generated by the operators D = Sujduj — 2ifd^ and II = w^c^o; — ^ujLpd^, 
i.e. it is two-dimensional. It is enough for equation (Q to be integrable with the Lie method. 

The invariance algebra Ai induces only the subalgebra spanned the operator D. Therefore, 
n is a pure hidden symmetry operator of the initial equation Let us note that the first 
non-trivial example of hidden symmetries connected with reduction of PDEs was found by 
Kapitansky ^1 for the Navier-Stokes equations. Wide classes of hidden symmetries of the 
Navier-Stokes equations were constructed in j^. See also for different notions of hidden 
symmetries of ODEs and other references therein. 

To reduce equation to an integrable form, we change the unknown function as 

, . 1 

ip = , I.e. ip = ^-^. 

In view of (@J) , V'tj 7^ and V' satisfies the equation iuj'^il^u?)^) = (V" "'^)a; which is simply integrated 
once. Further integration of the obtained first-order ODE oj'^iI^J^ = + Ci, where Ci is an 
arbitrary constant, with separation of variables results in the implicit solution 
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I , = — ± 1- Oq. 
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The integration constant Ci is normalized to {—1,0,1} by induced scale transformations. Let 
us note also that all values of Cq are equivalent with respect to the hidden symmetry group 
generated by the operator 11. Calculation of the integral depends on the value of Ci. 

If Ci = then ip > 0, i.e. in view of the definition of solutions may exist only for negative 
values of u. Integration under the condition Ci = results in the solution ip = (Co — |x~^)3 
corresponding to solution 1) of equation The integration constant Co is normalized to e G 
{—1,0, 1} by scale transformations associated with the operator D. 

The condition Ci = — 1 implies < ■0 < 1, i.e. solutions again exist only for negative values 
of uj and have an implicit form giving solution 6) of equation . In case Ci = 1 we have the 
constrain 0<— lor'(/'>0 and derive solution 7). 

The other kind of Lie invariant solutions are the similarity solutions which are invariant with 
respect to scale transformations. It is more convenient here to work in terms of the variables 
(t, X, v). The ansatz constructed with the subalgebra (D) has the form v = tip{uj), where uj = x/t. 
After substituting it to (jHl), we obtain the reduced ODE 

f^LULu - T;{Vi^) + ^fu; - (p = 0. 

6 2 UJ 

It has two polynomial solutions ip = —\uj and if = |a;^ — |a; which correspond to solutions l)e=o 
and 2) of equation ©. 



5 Non-Lie ansatz 

The form of Lie invariant solutions 1) and 2) directs us to look for more general polynomial 
solutions of equation (jSJ. As a result, we find the ansatz 

V = Lp^{t)x^ + (p^(t)x^ - ^x, 
which reduces equation Q to the system of two ODEs 

Up to translations with respect to t and scale transformations induced by the Lie symmetry 
group of the initial equation, the above system has the following inequivalent solutions 

(e^e), (0, -ft), (l,-3tan2t), (-1, 3tanh 2f), (-l,3coth2i) 

which correspond to solutions 1)^5) of equation Solutions 3)-5) are non-Lie ones. 
The above ansatz is rewritten in terms of the function u as 

/ 9 \ -2/3 

u = [if^ {t)x^ + (f'^ {t)x^ ~ 4 ^ ) 

This ansatz can be interpreted in the framework of a number of different approaches of find- 
ing exact solutions of nonlinear PDEs, such as nonlinear variable separation P|, the method of 
differential constraints (22|) anti-reduction [HI or generalized conditional symmetries I27j . 
(See also |16j for connections between these approaches.) Thus, the differential constraint 
2x^{x~'^u~'^/'^)xx = —9 corresponding to the ansatz is compatible (i.e. in involution) with equa- 
tion ((2j. 'Anti-reduction' of equation ((2)) by the ansatz containing two new unknown functions 
of one argument to the system of two ODEs means that 

(8n^ + Sxuux + Bx'^Ux'^ — x^uu^x + &xu'^'^)du 
is a generalized conditional symmetry operator of equation (jSJ. 
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6 On nonclassical symmetries 



We also study nonclassical (conditional) symmetries of equation @. (See e.g. [2^1121] for neces- 
sary definitions and properties of nonclassical symmetries.) Reduction operators of equation @ 
have the general form Q = rdt + ^dx + ??5n, where r, ^ and r] are functions of t, x and u, and 
(r, ^) 7^ 0. Since ^ is an evolution equation, there are two principally different cases of finding 
Q: r / and T = 0. 

We derive the system of determining equations in case of r 7^ and integrate it completely. As 
a result, we obtain the following statement. Any conditional symmetry operator of equation (j2]) 
in case of t ^ is equivalent to a Lie symmetry operator. 

As well-known, the operators with the vanishing coefficient of dt form so-called 'no-go' case 
in study of conditional symmetries of an arbitrary (1 + l)-dimensional evolution equation since 
the problem on their finding is reduced to a single equation which is equivalent to the initial 
one (see e.g. j71|23HSj). Note that "no-go" has to be treated as impossibility only of exhaustive 
solving of the problem. A number of particular examples of reduction operators with r = can 
be constructed under additional constraints and then applied to finding exact solutions of the 
initial equation. Since the determining equation has more independent variables and, therefore, 
more freedom degrees, it is more convenient often to guess a simple solution or a simple ansatz 
for the determining equation, which can give a parametric set of complicated solutions of the 
initial equation. Namely, in the case r = we have ^ 7^ 0. Up to usual equivalence of reduction 
operators, ^ can be assumed equal to 1, i.e. Q = dx + rjdu- The conditional invariance criterion 
implies the determining equation on the coefficient t] 

_ r]xx + '^Vyxu + v'^Vun _ + 6t?^t?„ 15??^ T] - ur]u _ _u_ 

~ m3/2 2^5/2 + 4m7/2 ^ x 

which is reduced with a non-point transformation to equation where rj becomes a parameter. 
We have found partial solutions of the determining equation but all of them result to the above 
solutions of equation @. In particular, the operator 

dx - 2^ (1 ± VCiu - xuV^ ) du 

gives solutions 1), 6) and 7) in case of the values Ci = 0, Ci = —1 and Ci = 1 correspondingly. 
The ansatz of Section |S1 is associated under the condition ip^ = Ci with the operator 

dx-— (4Cix^ni + 9xu^ + 8)du 
6x 

which results, therefore, in the solutions l)-5) depending on values of Ci. 
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